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Abstract 

The adsorption kinetics of the cationic surfactant dodecyltrimethylammonium bromide at 
the air-water interface has been studied by the maximum bubble pressure method at concen- 
trations below the critical micellar concentration. At short times, the adsorption is diffusion- 
limited. At longer times, the surface tension shows an intermediate plateau and can no longer 
be accounted for by a diffusion limited process. Instead, adsorption appears kinetically con- 
trolled and slowed down by an adsorption barrier. A Poisson-Boltzmann theory for the elec- 
trostatic repulsion from the surface does not fully account for the observed potential barrier 
The possibility of a surface phase transition is expected from the fitted isotherms but has not 
been observed by Brewster angle microscopy. 
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Introduction 



Over the last decades, the adsorption kinetics of surfactants at fluid interfaces has been the subject 
of many studies due to its prime importance in numerous applications such as wetting, detergency, 
foaming and emulsification.- Adsorption kinetics is commonly studied by creating a freshly ex- 
posed surface in contact with a bulk surfactant solution and measuring the temporal dependence of 
its surface tension. This dynamic surface tension^-^ has been shown to correlate much better than 
the equilibrium surface tension with properties crucial for applications, such as foaming ability of 
surfactant solutions and the spreading velocity of fluid films on top of solid substrates.-^ 

Experimental work on dynamic surface tension stimulated over the years theoretical modeling. 
In their pioneering model from the 1940s, Ward and Tordai^ assumed that the adsorption is limited 
by diffusion (DLA), resulting in an asymptotic decay as the inverse square root of time. The 
Ward-Tordai model is still widely used nowadays to analyze dynamic surface tension data. Yet, 
the improved accuracy achieved in recent years has revealed important deviations from the DLA 
behavior in several systems.^ These deviations have been attributed to a number of factors such as 
the existence of kinetic adsorption barriers, especially for charged systems;^ the role of surface- 
active impurities;^ lateral relaxation modes in the monolayer plane; ^ and the formation of surface 
crystals.— We note that the effect of electrostatic barriers on the kinetics of surfactant aggregation 
was analyzed also in the different context of ionic micelles.— 

For example, an anomalous long time decay with an intermediate plateau has been observed 
for the ionic solutions of sodium dodecylsulfate (SDS) with no added salt. Since surfactant hy- 
drolysis produces trace amounts of dodecanol that slowly adsorb at the water-air interface, some 
authors attributed the peculiar surface tension kinetics to these surface-active impurities.^ A simi- 
lar anomalous behavior was reported for aqueous solutions of aerosol OT,^'^^ as well as for SDS 
at the alkane-water interface. In the latter case, dodecanol impurities are not expected to accu- 
mulate at the interface and to affect the adsorption dynamics because they are soluble in the alkane 
phase.— 

While the origin of non-DLA-like behavior is not well understood for all systems, there is 
strong evidence that the adsorption of pure nonionic surfactants is usually diffusion-limited, 
For ionic surfactant solutions, however, an electrostatic surface potential is progressively created 
as the adsorption proceeds. This surface potential acts as an adsorption barrier for additional 
surfactant molecules as they migrate from the bulk toward the surface, thus potentially giving rise 
to a non-DLA behavior. The dynamic surface tension of ionic surfactants (SDS) at the oil-water 
interface was previously shown to be different from DLA behavior,'^ but when salt was added to 
the ionic surfactant solution, the process returned to a DLA-like behavior. This can be explained 
by electrostatic screening of the surface potential in the presence of added electrolyte. 

In the present work we extend the results of Ref [l7| to the adsorption of ionic surfactants at the 
air-water interface. We have studied aqueous solutions of the cationic surfactant dodecyltrimethy- 
lammonium bromide (DTAB) at various concentrations below the critical micellar concentration 
(cmc) using the maximal bubble pressure (MBP) method. At short times the kinetic behavior is 
found to be DLA-like, while at longer times the kinetics appears to be slowed down by an adsorp- 
tion barrier. The experimental results are analyzed and compared with the theoretical predictions 
of Ref 15. The influence of the adsorption barrier is better seen just before the critical micellar 
concentration (cmc) and working below the cmc has the advantage that the theoretical interpreta- 
tion is more straightforward as it does not have to take into account the presence of micelles. In the 
SDS case, hydrolysis produces dodecanol, which is very surface active at these concentrations (at 
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oil-water interfaces it is less a problem, because dodecanol is solubilised in the oil phase). We have 
chosen here to work with DTAB instead of SDS, because DTAB is much more stable chemically. 



Materials and Methods 

Dodecyltrimethylammonium bromide (DTAB) was supplied by Sigma (> 99% purity) and used 
without further purification. De-ionized water was obtained from a DEIONEX MS 160 equipment, 
with a resistivity of 12 MQ. All the measurements were performed at 25 ±0.2°C, and the solu- 
tions were prepared by dilution of concentrated stock solutions. The equilibrium surface tension 
measurements were conducted using a circular Teflon trough (capacity 5 ml), housed in a Plexi- 
glas box with an open-frame version of the Wilhelmy plate. We used a rectangular open frame 
(20 X 10 mm), made of a platinum wire and attached to a force transducer mounted on a motor, 
allowing it to be drawn away from the surface at a controlled constant rate. 

The dynamic surface tension was measured bya maximum bubble pressure (MBP) method. 
This in-house-made instrument (described in Ref |20|) is used to measure the maximum pressure 
necessary to detach a bubble from a capillary. The surface tension y is then obtained from the 
Young-Laplace equation, y = RAP/2, where AP is the maximum pressure difference between the 
gas in the bubble and the surrounding liquid, and R is the internal radius of the capillary. We have 
used a glass capillary with an inner diameter of 200 /im, hydrophobized with hexamethyldisilazane 
(Sigma). The advantage of such a device over methods based on analysis of the bubble (or drop) 
shape, is that it allows access to short adsorption times, down to about 10 ms, whereas the shape 
methods are limited to times longer than 1 or 2 s. The MBP instrument was tested with the non- 
ionic surfactant Triton XI 00, and the results were found to be compatible with DLA kinetics, 
indicating that uncontrolled convection effects were negligible even at times down to 0.1 s. 

With the MBP instrument, the bubble surface age is usually associated with the time interval 
between consecutive detachments of two bubbles. However, there is a dead time (delay) between 
the detachment of a bubble and the formation of the next bubble. In a recent paper, Christov et 
al.— showed that the effective aging time of the interface, t, is smaller than and proportional to 
the time actually measured between consecutive detachments, ?age: t = tage^^^, where A is an 
apparatus constant described in more details in Ref From the experiments with Triton X-100 
and its diffusion constant (D = 4 x lO^^cm^/s) the apparatus constant, A = 4.3 ±0.3, is obtained. 
In the following, all the curves are represented as a function of the corrected time, t. 

Since it is known that different MBP experimental setups give different dynamic surface tension 
curves 7(?)>~ experiments were also performed with a commercial MBP instrument (MPT2 from 
Lauda) for comparison with the in-house made apparatus. We also employed a Brewster angle 
microscope (Mini-BAM from NFT) together with a Langmuir trough (NIMA 601 BAM) to image 
the surface monolayer and visualize any of its possible surface heterogeneities. 



Results 

[figure][l][]l shows the reduction in the equilibrium surface tension, A/eq = 7eq — 7w, for salt- 



free DTAB solutions as a function of the DTAB concentration c, where 7w is the bare air-water 
surface tension. The data obtained by the two methods — the Wilhelmy plate and the MBP at long 
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times (all measurements reach the equilibrium values for fage < 50 s) — are in good agreement. 
For concentrations above the cmc (= 15 mM), the surface tension saturates to a constant value, 
A/eq ~ -33 mN/m. 



[figure] [2] [] 2 shows dynamic surface tension curves for DTAB concentrations below the cmc. 



At the higher concentrations (panels d-f) the curves manifest an intermediate plateau indicating 
a double-relaxation process. The behavior observed for high DTAB concentrations at the air- 
water interface is qualitatively similar to that found for SDS at oil-water interfaces at lower SDS 
concentrations (the cmc is smaller for SDS). In the latter case, however, the characteristic times 
were longer (because the concentrations were smaller) and accessible by the drop shape method. 




c (mM) 

Figure 1: Equilibrium surface tensions obtained by tine Willnelmy plate (circles) and MBP (up- 
triangles) techniques. The solid line is a fit obtained by using eqQ] The values of the fit parameters 
are: a = 11.8, /3 = 7.6, and a = 7.2 A. The upper and lower dashed lines show the theoretical 
curves using the same values for a and a but with /3 = 7.2 and 8.0, respectively. 



The plateaus seen could correspond to the occurrence of a phase transition in the surface layer 
during which surface pressure n ~ y— remains constant.^-' In order to check for this possibility, 
we imaged, using a BAM apparatus, the surface at surfactant concentrations close (below and 
above) the "knee" (c ~ 2mM) in |[figure][l][]l[ but no monolayer heterogeneities were observed 
(images not included). Note, however, that this does not exclude the presence of domains with 
sizes below the optical resolution (f^ 5 x 10^^ m). 



[figure] [3] []3 shows a comparison of the dynamic surface tension of solutions containing 



10.8 mM DTAB without salt (same as [figure] [2] []2p and in the presence of 5 mM NaBr. Placing 
the two curves on the same figure shows that the small amount of added salt has almost sup- 
pressed the plateau. In fact, the line in that figure represents a fitting to a DLA process using 



eq [[equation] [2] [] 2 
D 



From this fitting and eq [equation] [4] [] 4 we obtain a diffusion coefficient 
7 X 10^^ cm^/s, which is rather close to the previously measured one,^^ D ^ 6 x 10^^ cm^/s. 



The difference between these two values of D is compatible with uncertainties in the apparatus 
constant, X . 



In [figure] [4] []4 we show the dynamic surface tension for c = 3.8 mM obtained with the com- 
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Figure 2: Dynamic surface tension for DTAB concentrations of (a) c = 1 .08 mlVI, (b) c = 2.16 mlVI, 
(c) c = 4.32 mlVI, (d) c = 8.65 mlVI, (e) c = 10.8 mlVI , and (f) c = 13 mlVI. Tine lines are fits to 
eq|[equation][2][]2[ wfiile in tine insets tliat data points are fitted witin eq|[equation][5][]5 



mercial MBP apparatus (MPT2). The results are in qualitative agreement with our in-house appa- 
ratus. Note that the plateau can equally be observed, but the characteristic times are much shorter 
(as compared with | [figure] [2] [] 2| ) . For this commercial MBP the apparatus constant A seems to 
be lower than 1 (probably due to some uncontrolled convection). Therefore, we did not analyze 
quantitatively these results. 
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Figure 3: Dynamic surface tension for a mixed solution of 10.8mlVI DTAB and 5mlVI NaBr (open 
circles), compared with the corresponding salt-free DTAB solution (open squares). The solid line is 
a fitting curve using the DLA model of eq|[equation][2][]2 



Data Analysis 



Equilibrium Data 

The equilibrium adsorption isotherm and equation of state for ionic surfactant solutions were pre- 
viously derived.-i^*^ They are given by 



i,eq 



Are 



eq 



kBT 



ln(l-<^)o,eq) + f0a 



2 T-u,eq 



ib(l>o.. 



eq) 



+ 1-1 



(1) 



In eq[il ^ = a^c is the bulk surfactant volume fraction (a being the average size of a surfac- 
tant molecule), ^o.eq is the surfactant area fraction (surface coverage) at equilibrium, and k^T the 
thermal energy. The parameters a and /3 are the Langmuir adsorption parameter and the Frumkin 
lateral interaction parameter, respectively (both given in units of k^T). Finally, the parameter 
b{a,^) = [;rZB/(2a^)]^/2, where the Bjerrum length, Zb = e^/{£kBT) ~ 7 A, characterizes the 
strength of electrostatic interactions, and £ is the water dielectric constant. 

We numerically solve the isotherm of eq[T]for ^o,eq and substitute the result in the equation of 
state to calculate Ajeq. |[figure][l][]l| shows the fit of A7eq(c), thus obtained, to the equilibrium 
data (up to the cmc). There are three parameters in the fit: a, /3 and a. The theoretical curve shown 
[figure] [l][]T]corresponds to a fitting with a = 1 1.8 ±0.2, j6 = 7.6 ±0.4, and a = 7.2 ±0.2 A. In 



in 



[figure] [1][]1 we also demonstrate the sensitivity of the fit to the value of the Frumkin lateral inter- 
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Figure 4: Dynamic surface tension for 3.8mlVI DTAB solution measured witin a commercial MBP 



apparatus (MPT2). The solid line is a fitting with eq |[equation][2][]2[ In the inset, the fitting is done 
with eq|[equation][5][]5[ Note that the equilibration times are about 10 times shorter than those of 



[figure][2][]2:. 



action parameter, /3 . The key point is that a relatively large value of /3 is required to reproduce the 
"knee" in A7eq(c) at intermediate concentrations. The need for such a strong surfactant-surfactant 
attraction at the interface for ionic surfactants was recognized before, where it was suggested that 
it may be related to the adsorption of a small amount of counterions, which reduces the electrostatic 
repulsion, permitting the interaction of the non-polar parts of the surfactant chains. 

When /3 > 7.5 the fit of the model with the equilibrium data predicts a surface phase transition 
above a certain critical surface coverage. Since the surfactant molecules are soluble in the bulk 
solution, the surfactants at the surface can be treated as in a grand-canonical ensemble, and at 
equilibrium there cannot be a coexistence region between dilute and dense domains. At a certain 
value of the chemical potential (or equivalently of 0b = c?c) 0o,eq should jump discontinuously 
without a change in 7. Kinetically, however, this increase in the coverage should proceed via 
nucleation and growth of domains of the denser phase. For the above-mentioned values of param- 
eters a, /3 and a, the transition is between 0o — 0-52 (dilute) and 0o — 0-74 (dense) and occurs 
at concentration c ~ 2.27 mM, which is consistent with the "knee" observed in the equilibrium 
surface tension ( [figure] [1][]1 1. This is also close to the concentration at which the dynamic sur- 
face tension curves deviate more clearly from DLA behavior, [[figure] [5][]5| shows the theoretical 
equilibrium coverage as a function of concentration, as calculated by numerically solving eqiUfor 
the same parameter values used in the fit of [figure][l][]l The range of surface coverage and bulk 
concentrations corresponding to the transition is rather small, and, furthermore, sensitive to the 
fitted value of j8. (See dashed lines in [figure][l][]l ) This can possibly explain why a transition 
has not been detected in the BAM experiments. 
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Figure 5: Theoretically calculated equilibrium surface coverage, ^o,eq> as a function of bulk concen- 
tration, c. The curve is obtained from eq[i|with the parameter values as used to fit |[figure][1][]1 
The dashed line section indicates a region of discontinuous (first-order) phase transition for 0.52 < 

(|)0,eq < 0.74. 



Kinetics at short times 

At short times (?age < 5 s) but yet larger than T^, and for low surfactant concentrations, as is demon- 
strated in [figure] [2] []2[ the dynamic surface tension curves fit quite well the asymptotic time 
dependence of a DLA process.—*^ 



AY{t > Td) - Area 1 




(2) 



The fitted values for and A/eq are given in | [table] [!][]! 



Table 1: Fitted values for the diffusion time and the equilibrium reduction in surface tension Ajeq 



using eq [equation][2][]2 We have also included the value obtained from [figure][4][]4 Note that in 



this case the adsorption times were not corrected by the apparatus function A. 



c (mM) 


Td (ms) 


A/eq (mN/m) 


D X 10^ cm^ s^i 


2.16 


1.05 


-5.56 


2.3 


3.24 


0.36 


-9.68 


3.1 


4.32 


0.52 


-13.1 


1.2 


3.8 (MPT2) 


0.053 


-12.0 
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According to the theory— Td, which characterizes the relaxation of Ay in a DLA process, is 
given by 
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eq 



TtD 



(3) 



where D is the bulk diffusion coefficient of the surfactant molecule. From the known value of 
the surfactant concentration, its molar weight 308 g/1 and density 0.684g/l, we can estimate the 
volume fraction, 0b- At sufficiently high concentration, after the diffusion step, the surface tension 
reaches values corresponding to an almost saturated monolayer (constant slope region of the ex- 



perimental A/eq as a function of logc in [figure] [1][] 1 ). Therefore, for c > 2 mM we will assume 
that (^^, 



eq 



1 . This leads to 



D 



1 



(4) 



The value of the molecular size, a ^ 7.2 A was obtained in the previous section from the fit to 
the equilibrium isotherm and is about equal to the value estimated from the molecular volume, a^, 
where a^%h. The calculated D values are listed in |[table][l][]l| for c values between 2 and 4 mM. 
These D values are a bit lower than the known value for DTAB diffusivity in solution, D ^ 6 x 
10^^ cm^/s. This may indicate that the diffusion into the sub-surface region is slightly slower than 



the diffusion in the bulk. All the curves in [figure] [2] []2 clearly show that the adsorption kinetics is 
not controlled by diffusion at times larger than ? > 1 s as will be discussed in the following section. 



Kinetics at longer times 



As is evident from [figure] [2] [] 2 the dynamic surface tension curves at longer times deviate sub- 
stantially from the DLA behavior and cannot be well fitted with an inverse-square-root temporal de- 



cay. This becomes more pronounced as the concentration increases. In the insets of [figure] [2] [] 2 
we re-plot the data on a semi-logarithmic scale, demonstrating that the final relaxation to equilib- 
rium is exponential. 



7(0 - Teq ~ e 



(5) 



This relaxation is consistent with a kinetically limited adsorption (KLA), where the process is 



hindered by adsorption barriers, [[table] [2][]2| lists the fitted values of the relaxation time T^. These 



values are 2-3 orders of magnitude larger than the diffusion times listed in [[table] [!][]! 
From the theory-i^*^ we expect 

Tk = Tdcxp ( - a - /3 (^)o,eq + 2e vz/^B 7^) , (6) 

where Td has been defined in eq [equation] [3] []3[ and ^ — (i/zq + i//£,)/2 is the average of the 
equilibrium electrostatic potentials at the surface and sub-surface layers. From eqs [[equation] [3] []? 
and [ [equation] [6] [] 6] with 0o,eq — 1 we get e^ir/k^T = [a + /3 +\n{Ka'^Dc^x^)]/2. Using the fitted 
values of Tk ( l[table][2][]2| ), D = 6 x 10"^ cm^/s, and the fitted equilibrium values a2^11.8,/3~7.6 
and a ^ 7.2A, we obtain for the four higher bulk concentrations, respectively, four rather similar 
values for the average surface potential: exjf/k^T ~ 12.35, 11.73, 12.88, and 14.96. 
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Table 2: Fitted values for the kinetic relaxation time T^. 



c (mM) 


Tk(s) 


1.08 


1.06 


z. Id 


U.OD 


3.24 


0.23 


4.32 


1.24 


6.48 


0.27 


8.64 


0.05 


10.8 


0.33 


13 


1.3 



On the other hand, the Poisson-Boltzmann theory— yields 



(2^<l>0,eq)^ 



(7) 



where b{a, ^) has been defined below eq\T\ Taking ^o,eq — 1 and a ~ 7.2 A, we get for the four 
higher bulk concentrations, respectively, exifo/k^T ~ 8.4, 8.1, 7.9, and 7. Since the potential 
at the subsurface layer must be smaller than the surface one, < i/Zq, we should expect to find 
i//o > V- The values found for Xj/Q, based on the equilibrium Poisson-Boltzmann (PB) theory, are 
comparable, but smaller than the aforementioned potential barriers \j/ inferred from the dynamic 
surface tension measurements. This may reflect some inadequacies of the PB theory to account for 
all the experimental results as reported here. 



Conclusions 

In the present work we have analyzed the equilibrium and dynamic surface tension of DTAB at the 
water-air interface for several concentrations, all below the critical micellar concentration (cmc). 

At short times DTAB adsorbs in a diffusion-limited process (DLA) with a ~ temporal 
relaxation and a diffusion coefficient D 2^ 10^^ cm^/s. At longer times, the DLA behavior is fol- 
lowed by a kinetically limited adsorption (KLA) with an exponential relaxation, as predicted by 
the theory of Ref Ts" for salt-free surfactant solutions. For the higher concentrations, the dynamic 
surface tension exhibits an intermediate plateau, followed by a final, exponential relaxation occur- 
ring over time scales of several seconds. This behavior is similar to the one reported earlier for the 
adsorption of SDS at a water-oil interface in the absence of salt.— 

The experiments indicate that the adsorption of DTAB at the water-air interface undergoes a 
qualitative change in behavior in the concentration range, 2 < c < 3 mM. We could not detect 
the predicted phase transition using Brewster angle microscopy (BAM), but this could arise either 
because the surface domains are smaller than the optical resolution or because the transition occurs 
at concentrations other than those investigated (the predicted concentration range is very narrow. 



see [figure] [4] []4 1. 
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The kinetically limited relaxation is related to an electrostatic barrier created by the charged 
surface. The relevance of electrostatics is supported by the observed strong effect of added salt on 
the adsorption kinetics. Our analysis indicates, however, that the Poisson-Boltzmann theory does 
not fully account for the observed potential barriers. The discrepancy might originate from effects 
related to the finite size of ions concentrated close to the surface.— We note that such effects were 
included in previous studies,'^^*^^ where equilibrium measurements were well fitted by a van der 
Waals isotherm with a Stem layer of bound counterions.^^ Similar modifications of the current 
theory may also affect predictions concerning the occurrence of a surface phase transition during 
the adsorption process. 

Another mechanism that might be invoked in principle involves the possible formation of pre- 
micellar aggregates at concentrations closer to (but below) the bulk solution cmc. In a previous 
theoretical study, it has been shown that the adsorption kinetics in the presence of micelles ex- 
hibits an exponential relaxation, with a relaxation time that is related to the exchange of surfactants 
between aggregates and the surface. However, the exchange time for small chain surfactants is of 
the order of a microsecond,-2^i^ and is too fast to be seen in our experiments. Furthermore, one 
expects added salt to promote aggregate stability and, thus, to push the system toward KLA, in 
disagreement with our measurements. 

In conclusion, the current work brings further insight into possible mechanisms of ionic surfac- 
tant adsorption at fluid interfaces. Further theoretical and experimental work is necessary to fully 
understand this important phenomenon that controls the dynamic behavior of the interfaces. 
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